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We find that certain higher-order mappings arise as reductions of the integrable discrete A-type KP 
(AKP) and B-type KP (BKP) equations. We find conservation laws for the AKP and BKP equations, 
then we use these conservation laws to derive integrals of the associated reduced maps. 
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The search for discrete integrable systems has received considerable attention in the past decade. 
This has resulted in the discovery of integrable mappings of the second order, e.g., the Quispel- 
Roberts-Thompson (QRT) mapping, 1 and discrete Painleve equations. 2 Apart from second-order in- 
tegrable mappings, the results for higher-order integrable mappings are few. 311 Discrete integrable 
systems have applications to various areas of physics, such as statistical mechanics, quantum gravity, 
and discrete analogues of integrable systems in classical mechanics and solid-state physics. Here, we 
study a novel class of higher-order integrable mappings that have bilinear forms. 

As an example, we discuss the following 6th-order mapping: 

Dx n+ ?,X n+ 2X n+ yX n X n _yX n _2Xn-?> + ^n+2 x n+l x n x n—\ x n--2 

+Bx n+ iX n X n -i+C = 0. (1) 

(Here and below A,B,C,D are arbitrary parameters). How can we obtain integrals for this mapping? In 
ref. 3, a method for the construction of integrals was proposed and integrable third-order mappings that 
possess two integrals were obtained. However, this method is not applicable to higher-order mappings 
because this method uses some ansatz at first and requires the use of high-performance computers. 
If we consider 6th-order, 8th-order and higher-order mappings, this method does not work, because 
current computer power is not sufficient. 12 

In this Letter, we propose a systematic method of constructing integrals for a class of higher-order 
integrable mappings that can be transformed into a single bilinear form without the use of computers. 
Our proposed method is based on discrete bilinear forms related to the A-type KP (AKP) and B-type 
KP (BKP) soliton equations. Conservation laws for integrable partial difference equations have been 
studied in refs. 13 and 14. 
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Before we discuss conservation laws for discrete systems, let us briefly recall conservation laws for 
continuous systems. 15 To this end, consider a (scalar) partial differential equation (PDE) A[jc,kW] = 0. 
The conservation law of such a PDE is a divergence expression 

; ax J 

which vanishes for all solutions of the given system. It follows that there exists a function A (called the 
characteristic of the given conservation law) such that Y,j = AA. Similarly, the conservation law 
of a scalar partial difference equation A[n, u n ] = is the expression ~ id)Pj = 0, which vanishes 
for all solutions of the discrete system. (Here Sj is a unit shift in the tij direction, and A[n,u n ] denotes 
a smooth function depending on n, u n and finitely many iterates of u n ). Again there exists a function 
A such that 

Y,(Sj-id)Pj = AA. (2) 

i 

We will call A the characteristic of the discrete conservation law. 

Here, we give a list of characteristics of the discrete AKP and BKP equations. 
Discrete BKP equation 

The discrete BKP (Miwa) equation 16 is given by 

Alk+ 1 ,l,m%l+ 1 ,m+ 1 + B%k,l+ 1 ,m^k+ l,l,m+\ 

+Clk,l,m+\1k+l,l+l,m +D%l,m^k+lJ+l,m+l =0. (3) 

We have found the following 6 explicit rational characteristics for the discrete BKP equation: 

-l,/+l,m+l .m 

tkj,m+ 1 1jt,Z+ 1 ,m+ 1 tk,l+ 1 ,m tk+ 1 ,l,m^k+ 1 ,/+ 1 ,m^k+ 1 ,l,m+ 1 

,m+ 1 



Ai=, 

+D 
A 2 = B 

+D 
A 3 =C 

+D 
A 4 = C 



1k,l,m1k,l+ 1 ,m^k,l,m+ 1 ^k+ 1,1+1 ,m+ 1 ^k+ 1,1+1 ,m^k+ 1 ,l,m+ 1 
f2{l)lk+l,l-l,m+\ fl(t+ l)tfe,/+2,m 



1k+l,l,mtk+ 1 ,l,m+ 1 1k,l,m+ 1 ^k,l+ 1 ,m^k,l+ 1 ,m+ 1 ^k+ 1 ,/+ 1 ,, 
f2{l)lk,l-\,m h{l + l)tfe+l,/+2,»i+l 



tk,l,mtk,l,m+ 1 ^k+ 1 ,/,m ,m+ 1 1 ,m+ 1 ^k+ 1 ./+ 1 ,m 

h{rn)Xk+\,i+\, m -\ h{m + l)%u,m+2 



Xjfc+1 ,1+ 1 ,mtk,l+ 1 ,mTjt+ 1 ,l,m ^k,l,m+ 1 1 ,m+ 1 1 ,/,m+ 1 
fe{m)tk,ljn-\ fo{™+\)lk+\,l+l,m+2 



1k+ 1 ,l,m^k,l+ 1 ,m^k,l,m ^k+ 1 ,/+ 1 ,m+ 1 1 1 tk+ 1 ,/,m+ 1 
/4(0 X W-l,w+l /4(^ + l)^+l./+2.m 



1 ,/,m+ 1 ^k,l,m+ l^k.l.m ^k+ 1,1+1 ,m^k,l+ 1 ,m^i+ 1 ./+ 1 ,m+ 1 
^t(/)t*+l,/-l,m + l)t*,/+2,ro+l 



As 



1 ,l,m^k,l,m^k+ 1 ,/,m+ 1 x t,/+l ,m+ l^M-l ,m^+ 1 ,/+ 1 ,m+ 1 
f5{m)lk+\,l,m-l fs(m+ l)li,/+l,m+2 



X/t+l ,/+ 1 ,mXA:+ l,l,m^k,l,m ^k,l+l ,m+ 1 ^k,l,m+ 1 fjfc+ 1 ,/+ 1 ,m+ 1 
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-B 



f5{m)% k . 



l+l,m-l 



A 6 = B 



+C 



^k,[+l,m^kj,m^k+l,l+l,n 

-l,/+l,m 
tk,l+ 1 ,m+ 1 tk,l+ 1 ,mXk,l,m 

,m+l 



_ fs(™+ l)^+l,/,m+2 

1 ,/,m+ 1 tk,l,m+ 1 1£+ 1 ,/+ 1 ,m+ 1 

/ 6 (*+ 1)^+2,1 
1 ,l,m+ 1 1 J, m Xk+ 1,1+1 ,m+ 1 

/6(^+ l)X/t+2,/+l,m 



^k,l,m+l^k,l,m^k,l+l,m+l ^k+l,l+l,m^k+l,l,m^k+l,l+l,m+l , 

where fi(x),--- ,fe(x) are arbitrary functions. (If we choose fj(x) as constants or (— l)*, we can pro- 
duce 12 characteristics. Throughout this Letter, we only consider reductions in which all independent 
integrating factors are obtained taking fj (x) =constant.) From these characteristics, we can obtain the 
associated conservation laws, using eq. ((U). For example, Pi,Pz and P3 associated to Ai are 



fi(k) 



^2 ^k— 1,1+1 ,m+l tk+ 1 ,l,m 
tk+ 1 ,1+ 1 ,m ^k+ 1 ,l,m+ 1 



D 



2 tk— l,l,m^k+ 1,1+1 ,m+ 1 



^k,l+l,m^k,l,m+l 



AB 



-BD 

-AD 
P 2 = Mk)[AC 
P 3 =Mk)(AB 



tk- 1 ,1+ 1 ,m^k+ 1 ,/+ 1 ,m+ 1 
tk,l+l,m+ltkj+l,m 

X k - 1 ,l,mXk,l+ 1 ,m+l1k+ 1 ,/,; 
X k ,l ,rn^k,l+\,m^k,l,m+l 
1k+l,l,mtk-l,l,m+l 

1k,l,m+l*k,l .111 
1k+\,l,m1k-\,l+\,m 



AC 



1k-l,l,m+l1k+l,l,i 



CD 



1k-l,l+l,mtk+l,l,i 
1k,l,m1k,l+l ,m 



Ik- 1 ,l,m+ 1 ^k+ 1,1+1 ,m+ 1 



■ + 



BD 
CD 



^k,l,m^k,l,m+ 1 ^k,l+ 1 ,m+ 1 ^k,l,m+ 1 

1k,l,m1k- 1 ,m+l 1k+ 1,1+1 ,m+ 1 
^k,l+l,m+l^k,l+l,m^k,l,m+l 
^k-l,l,m^k+l,l,m+l ' 



.in 

^k-l,l,m^k+\,l+\„ 



"k,l ,mtk,l+l ,m ^k,l+l,m%k,l,m 

where P\, P2 and P3 correspond to the k-,l- and m-directions, respectively. 
Discrete AKP equation 

The discrete AKP (Hirota-Miwa) equation 16 ' 17 is given by 

Alk+lJ,mtkJ+l,m+l + BXkJ+l,mtk+l,l,m+l 
+Clk,l im +l1 k +l,l+l,m = 0. 



(4) 



Note that the discrete AKP equation is the special case of D = of the discrete BKP equation. The dis- 
crete AKP equation inherits the above 6 characteristics (with D = 0) from the discrete BKP equation, 
and we have found the following additional characteristic: 

fl{k+l+m)x k i m fi{k + l + m+\)x k+u+ i M+ i 
A 7 = , 

tk+ 1 ,l,m^kj+ 1 ,m"k,l,m+ 1 ^k,l+l,m+ 1 1k+ 1 ,/,m+ 1 1k+ 1,1+1 ,m 

where fj (x) is an arbitrary function. 

Reduction to Finite-Dimensional Mappings and Construction of their Integrals 
First example: Consider the following 4th-order mapping: 

Dx n j r 2X n j r j X n X n _\X n _2 + Ax n _|_ \X n X n _\ + B + C/x„ = 0. (5) 

Using the transformation of the dependent variable x n = x n+ \X n _i/x^ , we obtain the bilinear form 



Dx n +3^n-3 +AT„+2T„_2 +BX n+ iX n -i + Cxi = 0. 



(6) 
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This bilinear form is obtained from the discrete BKP equation by applying the reduction x„ = 
tzik+z 2 i+z 3 m, where Z\ = 1,Z2 = 2,Zj, = 3. Using the characteristics of the discrete BKP equation, 
we obtain the following integrating factors for the discrete bilinear form eq. (©: 

A B 

Ai = -A 6 , A 2 = -A4 , A 3 = -A4 + -A 6 , 
Aj= C|— X JI±l—) +A f T "- 4 %n+A 



^n+Y^n^n— 3 ^n^n— l^n+3 / 2^n— 3^n+l '^n+2'^n— l^n+3 

A 6 = gf T - 2 3*2 — Ucf Tii 1 



v ^«+2^n— l^n— 3 tn+lt/i— 2~n-\-3 / V^n^/i— 3^n+2 ^n^n— 2^n+3 , 

It is confirmed using the bilinear form eq. © that the integrating factors Ai, A2, A3 and A5 lead to 
the indicated linear combinations of A4 and Ag. Note that the two integrating factors A4 and A(, are 
independent. From the above integrating factors, we can write the integrating factors in terms of the 
;e- variable: 

A4 = T„_iT, 1+ iA 4 =c( ) +A (x„_ 2 *«-3 -X n+ T,X n+2 ) , 

A6 = T„_iT„ + iA6 

B\ ! +c 



, , 2 2 

Xn+lXnXn—lXn—2 *»+2*n+l*re*Jt— 1 / \ ^-n+l^n-^n— \%n— 2 ^n+2^ n j r \^n^n—\ 

We then obtain the following two integrals: 

Qa = CDx n+ 2X n+ lXfPC n — i ADx n j^T 1 X n j r 2X n j r iX n X n _iX n _2 

2, if 1 1 

—A \X n +T,X n ^2X n +{X n X n -\ + X n ^2X n +\XnX n -\Xn-2) — C I 1 

\Xn+2Xn+lXn ^n+l^n^n—l 
-BC + + - AB X n+ 3-jX n+2 -j —AC 2^ " " , 

( 1 111 

Q 6 = BDXn + 2X n+ lX n X n -\ +CD[x n+ 2X n+ lXn+X n+ \X n Xn-l) - AB I 1 1 1 

\%n+2 X-n+l %n %n—l 

/ 1 1 1 \ B 2 C 2 

AC + + 



Xn+2%n+\ %n+\Xn Xn%n—\ J Xn-\-2Xn-\-\Xn%n—l X n j r 2X n ^^X n X n —i 

-bc( L^ + . i ^ 



2 2 / ' 

\ Xn+2Xn J r\X r r^n— 1 X n ^2X n _^_iX n X n — i / 

It is not difficult to show that (?4 and <26 w& functionally independent. 

In the special case of D = 0, the fourth-order mapping eq. © reduces to the second-order mapping 

Ax n+ ix n x n -i + B + C/x n = , (7) 

which is a special case of the QRT mapping. 1 Using the transformation of the dependent variable 

Xfi — tn+ltn— 1 

/l n , we obtain the bilinear form 

AX„ +2 T„-2 + 5X„ + lT„-l + Cx 2 = . (8) 
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This bilinear form is obtained from the discrete AKP equation by applying the reduction x„ = 
tz l k+z 2 i+z 3 m, where Z\ = 1,Z2 = 2,Zt, = 3. Using the characteristics of the discrete AKP equation, 
we obtain the following integrating factors for the discrete bilinear form eq. ([8]>: 

Ai = , A 2 = -Ai , 

\'^n^n+2'^n— 1 ^n— 2^n^n+l / 

c c 2 

A 3 = -Ai , A 4 = BAi , A 5 = - —A, , 
C 

A 6 = -AAi , A 7 = -Ai . 

There is only one independent integrating factor, Ai. From Ai, we can write the integrating factor in 
terms of the x-variable: Ai = x n +iT„_iAi = 1 /x n+ \ — 1 /x n -\- We then obtain the following integral: 

( \ 1\ C 

Qi = -Ax n+l x n +B + — + 



K Xn+l Xn J X n -^iX n 

Second example: As a second example, let us discuss the following 6th-order mapping: 

Dx n+ 3X 2 + 2 x ^+l x n x n-\ x n-2 x n-3 + Ax n+2 x\ + yx\x\_\X n - 2 + Bx n+ iX n X n -i +C = 0. (9) 

Using the transformation of the dependent variable x n = x n +iX n _i/x„ , we obtain the bilinear form 

AT„+3T„_3 +Bx n+2 X„-2 + CX n+ iX„_i + Dx„ +4 x n _ 4 =0. (10) 

This bilinear form is obtained from the discrete BKP equation by applying the reduction x„ = 
tzik+z 2 i+z 3 m where Z\ = \,Z 2 = 2,Zt, = 5 or Z\ = l,Z 2 = 3,Zj = 4. Using the characteristics of the 
discrete BKP equation, we obtain the following integrating factors for the discrete bilinear form eq. 
CEO): 

Ai = -A 6 , A 2 = A4 , A 3 = -As , 
A 4 = C ( ^-1 ^ I + ,\ 





-1 






V 


hl 




~n+2tn-{ 


-lV 


-4 




-lV 


-2V 


h4 


X n - 


-8 








-8 




"n—Y^n- 


-3V 


-4 


**n- 




-lV 


^4 


tn- 


-3 






t/H 


-3 




»n+3*n- 


-2^n- 


-4 


»n- 


' r 2^n- 


-3V 


h4 



A 5=A 2-2 ^ ) +B 



A 6 = B[ — ^ )+C 













-3^n— 4^n- 


h2 ^n+3^n— 2^n- 


h4 










Xjj- 




-1 ^1+2^+1^ 


-4 




X>] 


X n 






hl^n-4^n- 


r3 ^fj— l^B— 3^fl- 


h4 



It is confirmed using the bilinear form eq. (iTOl) that integrating factors Ai, A 2 and A3 lead to A$, A4 
and A5 respectively. Note that the 3 integrating factors At, A5 and A6 are independent. From the above 
integrating factors, we can write the integrating factors in terms of the x-variable: 

A4 = x n+ iX, l _iA 4 

C ( 1 1 \ A f 



Xn+lX 2 X n —[ \x n _^X n _2X n —T, x n+3 x l l ^- 2 X n-\-i I x n+\ x n x n—\ 



A5 — x„+iX M -iA5 — Ax n (x n _YX n _ 2 x n _^x n _^x n _^x n _( y x n --i x n Yl x n+6 x n+5 x n+A x n+3 x n+2 x n+\) 



+Bx n (x n _ iX„_2 x n-3 x «-4 x n-5 x n-6 x n+6 x n +5 x n+4 x n+3 x n+2 x n 



+V> 



5E 



J. Phys. Soc. Jpn. LETTER 

A A B /ll 



2 2 2. 
x n+2 x rl J r \ x n X n-\ Xn -'2- \ x n— 3 x n+3 

1 1 

+ 



x n+2 x l + i x n x l-l x n-2 \ x n-\ x n-2 x n-3 x n+3 x n+2 x n+\ 

We then obtain the following three integrals. (Note that eq. (O can be used to eliminate, e.g., x n+ 4 and 
x„_ 3 from Q 4 , and similarly for Q 5 and Qe-)'- 

Q.4 = CD(x n+ ^X^ l+ 2 x n+l x n + x n+2 x n+l x n x n- 1 ~\~ x n+\ x ^i x ^i-\ x n-2) 

ADx n ^4X n _i_? l X n _i_2 x n +i x n x n — \ x n—2 x n— 3 ^ x n+3 x «+2 x rc+l x n x rc— l x n— 2{ x n+4 x n— 3) 

2 : 6 

— BC £ 2 2 ^ 12 x n+4-j x n+3-j 

j=0 X n+3-j x n +2-j X n+\-j X n-j j=0 

-C 2 ( 2 3 2 + 2-4l V AC I X " +4 -^ +3 ^' , 

yx n+ 3X n+2 X n+l X n X n -i x n+2 x n+ i x n X n-l X n-2 J j =0 x n+l-j x n-j x n-l-j 

3 

Q5 = AD £ x rc+7-j x; n+6-j x; n+5-j x; n+4-j x; n+3-j x; n+2-^ 
7=0 

2 

-\-BD £ X H+6-y x n+5-y x n+4-j' x «+3-j x «+2-j x «+l-j'^ 
7=0 

4 

+A £ X n +7_yX n+6 _yX n+5 _yX n+ 4_yX ;1+ 3_^X n+ 2-^„+l-^„-j- :li; , 1 -l-j- :li; «-2-y 
7=0 

4 

+AB £ x n+6-;^„+5-y x „+4-^«+3-y^«+2-7 X n+l-7 X fi-j X fi-l-y' x n-2-7 
7=0 

4 

XL X n+ 7 -7 a n+6-; X /)+5-7 X n+4-y x n+3-y x n+2-y x n+l-y' x «-y' x n-l-7 
7=0 

5 

+B £ *n+6— yX n+ 4_ yX n+ 3_ jX n+ 2—j x n+\—fhi— j x n— 1—j 
7=0 

5 

+SC £ X n+ (>-jX n+5 _jX n+4 _jX n+ T l _jX n+ 2-j x n+l-j x n-j 
7=0 

6 

+AC £ x n+7-j x n+6-j x n+5-j x n+4-j x n+3-j x n+2-j x n+\-j x n-j 1 
7=0 

3 

26 = BDx n+ T,X n+ 2 x n+\ x n x n-\ x n-2 + CD £ x n+3- j x n+2- j x n+\- j 

7=0 

5 1 2 1 

-AB £ -i AC £ l - 

y'=0 ■*n+3— 7 ;=o x «+3— 7' 1| -«+2— 1 —j x n—j 

B 2 C 2 



2 3 3 2 
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BC 



BC 



2 2 2 .2 2 2 

x n+3 x n+2 x n -^ \ x n X n — \ x n— 2 x n+'i x ^ l -L-'2 X n+ 1 x n x n— \ x n— 2 

Using a Mathematical software package such as Mathematica, one can show that Q4, Q$ and Q(, are 
functionally independent. 

We consider the 4th-order mapping 



2 2 2 

^ x n-\-2 x n +\ x n x n-\ x n-2 + Bx<n+\ x n x n- 1 +C = 0. 



(ii) 



This mapping is the special case of D = of the 6th-order mapping eq. (©. Applying x n 
/x„ , we obtain the bilinear form 



AX„ +3 X„_3 +BX„ +2 X„_2 +CX„ + \% n -\ = 0. 



(12) 



This bilinear form is obtained from the discrete AKP equation by applying the reduction x„ = 
tz [ k+z 2 i+z 3 m, where Z\ = 1,Z 2 = 2,Z^ = 5 or Z\ = 1,Z 2 = 3,Z3 = 4. Using the characteristics of 
the discrete AKP equation, we obtain the following integrating factors for the discrete bilinear form 
eq. (fT2l : 

X M +2 X„_2 



Ai 

A 2 



Xn X w 



X«— 3Xn+2X«+l X n — 2Xn+3X)!— I 



c z 



A 3 = "X2 A 1 

c 3 

As = ^ Ai + 



B 2 C 



A 3 
B 2 C 2 



-A 2 , A4 = -BA 2 , 
A 2 , A 6 = -AAi , 



A 7 



C 



-A 



2 • 



It is confirmed using the bilinear form eq. (fT2l that the two integrating factors Ai and A2 are in- 
dependent. From the above integrating factors, we can write the integrating factors in terms of the 
x-variable: 



Ai = x n+ iX„_iAi = 

A2 = X n+ iX„_iA2 = 

We then obtain the following two integrals: 

1 



1 



1 



1 



x n+\ x n x n— 1 \ x n+2 x n—2 



1 



1 



1 



x n+\ x n x n—\ \%H—l%n—2 x n+2 x n+\ 



Qi — —Ax n+ 2X n+ \x n x n _\ +B 



1 1 1 

+ + — + ■ 



Q2 — A { x n+2 x n+l + x n+l x n + x n x n- 1 , 



x n+2 x n+i x n x n—\ 
1 



+ ■ 



c 



x n+2 x n+\ x n x n—\ 



B 



+ ■ 



1 



\ X n+2 x n+l x n x n+\ x n x n—\ 

We note that eq. (fTTT > preserves the symplectic structure 



C 



2 2 

x n+2 x n -\- 1 x n x n— 1 



x n—2 x n—\ 
~ x n—\ x n—2 
x n x n—2 x n x n— 1 



Xn—$%n x n—2 x n+\ 
x n—\ x n x n—l x n+\ 
X n X n -\- 1 



V x n+\ x n—2 x n+\ x n—\ x n+\ x n 
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and that the two integrals Q\ and Q2 are in involution w.r.t. this structure, giving an independent 
confirmation of the integrability of eq. (flTT) . 

We have studied a class of integrable mappings that can be transformed into a single bilinear form. 
We have proposed a method of constructing integrals of these higher-order integrable maps. The key to 
the construction is the conservation laws of the discrete bilinear forms of the associated AKP and BKP 
equations. Note that, generalizing the examples in this Letter, we can construct a family of higher-order 
mappings from the discrete AKP and BKP equations, by applying the reduction x„ = lzik+z 2 i+z 3 m for 
any Z\, Z2 and Z3. 

In general, soliton equations in the 1+1-dimension have infinitely many conservation laws. One 
may expect the existence of infinitely many conservation laws for discrete AKP and BKP equations. 
However, finding infinitely many conservation laws in the form of eq. (2) is difficult in the case of 
2+1 -dimensional discrete soliton equations. Finding infinitely many conservation laws for discrete 
AKP and BKP equations is still an open problem. 

We hope to discuss details of our methods and higher-order mappings in the class given here in a 
forthcoming paper. 
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